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miTFJIODYNE OPTICAL SPECTRUM ANALYZER 

CLAIM OF BENEFIT OF PROVISIONAL PATENT 
APPLICATION 

Priority is claimed from provisional application 60/394,261 filed 
5 on July 9, 2002. The provisional application is incorporated by 
reference. 

RELATED APPLICATIONS 

This appUcation is related to conmonly-assigned application 
serial number 10/ , , entitled "Polarization Diversity Detection 
10 Witiiout a Polarizing Beam Splitter " filed on July 8 , 2003. 

FIELD OF THE INVENTION 

The present invention relates to optical measurements, and more 
particularly, to a method and system for optical spectrum analysis 
utilizing optical heterodyne detection. 

15 BACKGROUND AND SUMMARY OF THE INVENTION 

Optical heterodyne systems analyze the optical spectrum of an 
input optical signal under test. A basic optical heterodyne detection 
system includes an optical coupler that combines an input optical signal 
from a fiber or device under test with a reference optical signal. The 
20 resulting "mixed" optical signal includes a heterodyne "beat" signal at a 
frequency equal to the frequency difference between the input optical 
signal and the reference optical signal. The reference signal frequency 
is changed or "swepf ' across the bandwidth of the input signal. The 
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detected beat signals over that swept bandwidth are processed to 
determine one or more characteristics of the input optical signal such as 
frequency, wavelength, or amplitude. 

Typically, heterodyned spectrum analyzers detect amplitude 
5 spectra at high resolution compared with other forms of spectrum 

analyzers. This difference in resolution can be as large as four orders of 
magnitude. Unfortunately, the local tunable laser source for 
heterodyned spectrum analyzers is very expensive, and the need for very 
fine accuracy of the instruments is not clearly shown. As a result, there 
10 are very few heterodyned spectrum analyzers on the market. 

A problem with conventional optical spectrum analyzers is they 
only provide ampUtude information about the input optical signal's 
frequency spectrum— phase information is not obtained. But absent that 
phase information, the input optical signal cannot be accurately 
15 reconstructed in the time domain. Indeed, two different signals 'may 

have the same amplitude spectrum but different phases. In this situation, 
a detected amplitude spectrum, by itself, cannot be resolved accurately 
into one of those two signals. 

This difficulty is illustrated in Figs. 1-5. Fig. 1 shows a binary 
sequence corresponding to an original time varying signal. The vertical 
axis is labeled the modulating signal and is graphed against time on the 
horizontal axis. The original time varying signal modulates an optical 
"carrier" resulting in a modulated optical signal, which when transmitted 
over a fiber, may be detected as a frequency spectrum as shown in Fig. 
2. One might assume that the original signal could be simply recovered 
by transforming the frequency spectrum into the time demand, i.e., by 
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applying an inverse Fourier transfonn to the detected signal. But this 
assumption is not true. The detected spectrum includes only the 
amplitude of the signal; the phase is missing. The origmal time-domain 
signal cannot be recovered in the inverse Fourier transform operation 
5 because the phase information has been lost. The resulting time domain 
waveform in Fig. 3 resulting from the transform is distorted beyond 
recognition from what is shown in Fig. 1. Indeed, it is compressed into 
a single pulse rather than the series of pulses shown in Fig. 1. 

If both the phase information as well as the ampUtude 
10 information can be determined from the frequency spectrum shown in 
Fig. 2, both the real and imaginary parts of the spectrum could be 
determined as shown in Fig. 4. Transforming the frequency spectrum of 
the signal that contains both real and imaginary parts, (both amphtude 
and phase information), permits recovery of the original information as 
15 shown in Fig. 5 which corresponds to the modulating signal information 
shown in Fig. 1. 

An objective of the present invention is to detect an optical 
signal's complex power spectrum, such as that shown in Fig. 4, and to 
determine from that detected spectrum both amphmde and phase so that 
20 an original information signal can be reconstructed accurately in the 
time domain. A heterodyne optical signal analyzer in accordance with 
the present invention performs such a signal reconstruction to accurately 
recover the full, complex (amphmde and phase) time domain signal as 
well as the spectrum, thus providing an instrument with greater utihty. 

25 All laser transmitters introduce an unwanted phase shift on the 

optical signal. This unwanted phase shift, referred to as "chirp," Umits 
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the distance over which the optical signal can be sent. The present 
invention enables precise characterization of the signal chirp permitting 
accurate prediction of the transmission distance limit, and even how to 
modify the transmitted optical signal so that longer transmission 
5 distances can be achieved. 

In addition to resolving both the phase and amplitude of an 
optical signal as a function of time and chirp characterization, the 
performance of the heterodyne optical signal analyzer in accordance 
with the present invention is not limited by the bandwidth of its 

10 components. For example, a 40 GHz optical input signal can be 

accurately processed even though the heterodyne optical signal analyzer 
operates at a much lower frequency. The heterodyne optical signal 
analyzer slowly builds a picture of the input signal by sweeping the 
reference signal below and above the frequency of the input signal to be 

15 reconstructed. As a result, signals may be reconstructed with very high 
resolution. Effective sampling rates as high as 1000 GHz may be 
achieved. The measurements produced at these rates do not have any 
detector bandwidths incorporated into them, and thus the signals are 
nearly free of distortion. 

20 Accordingly, the heterodyne optical signal analyzer in 

accordance with the present invention permits optical communication 
system designers to determine in advance what an optical signal will 
look like when it is transmitted over a fiber network. That information 
is extremely important to people who design such networks. Currently, 

25 no instrument is available that permits one to identify the optical signal 
being propagated at any point in the network. In other words, the 
transmitted signal is an unknown. Although the signal's power as a 
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function of time and the amount of chirp (variation in frequency as the 
power is tumed on and off) can be measured, how the amplitude and 
phase of the signal vary with time is unavailable. Until now, there is no 
reliable way to predict how a particular transmitter will interact with a 

5 specific set of optical components or an overall optical link. Knowing 
the precise nature of the signals present in a network is the first step in 
making the network cheaper and more efficient. With conventional 
optical signal analysis instrumentation, there is a great deal of 
interpretation and guess-work involved leaving simple analysis 

10 problems to highly skilled and expensive people. The invention reduces 
the problem of analyzing an optical network to an analysis closer to that 
employed for existing cable or microwave networks. 

Other aspects of the present invention relate to caUbrating the 
heterodyne optical signal analyzer. Light in a standard optical fiber 

15 includes two orthogonal polarization modes. To fully characterize the 
light in the heterodyne optical signal analyzer, a vector representation of 
the light requires determining amplitude and phase for each of the two 
orthogonal modes. The components of a heterodyne optical signal 
analyzer, including optical couplers, detector blocks that detect optical 

20 power and convert it into an electrical signal, and the reference signal 
generator, all have errors and offsets. For example, optical power 
detectors are very sensitive to changes in polarization of the optical 
input signal and the reference optical signal. Several different 
calibration procedures, including detector calibration, vector calibration, 

25 and reference signal caUbration are described below. 

In a general example embodiment of the invention, an optical 
signal analyzer includes a first coupler, a first detector block, and a data 
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processor. The first coupler mixes an optical reference signal and an 
optical input signal whose characteristics are to be determined and 
generates multiple mixed signal outputs. The first detector block detects 
multiple power signals from the multiple mixed signals. Each 
5 individual detector in the detector block includes in one example 
implementation a photodetector, an amplifier, an analog-to-digital 
converter for converting the amplified output into digital power signal 
information, and a buffer for storing the digital power signal 
information. 

10 The data processor determines the original optical input signal in 

the time domain from those multiple detected power signals. Both the 
ampUtude and the phase of the optical input signal are determined from 
the detected multiple power signals for each different frequency of the 
reference signal as it is swept across the frequency bandwidth of the 

15 input signal. From these detected outputs over the swept frequency 
range, the original time domain signal is reconstructed using signal 
processing procedures outhned below. 

In the first example embodiment, the first coupler generates first, 
second, and third mixed signals, and the first detector block detects 

20 corresponding first, second, and third power signals. But in a second 
example embodiment, a second coupler and a second detector block are 
added to take into account the two polarizations of tight, thereby 
permitting more accurate signal reconstruction. The first coupler and 
detector detect mixed signals where the reference signal has a first 

25 polarization. The second coupler and detector detect mixed light where 
the reference signal has a second different polarization. The second 
coupler generates fourth, fifth and sixth mixed signals, and the second 
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detector detects foxirth, fifth and sixth power signals from the fourth, 
fifth and sixth mixed signals. 

In the second example embodiment, the data processor 
determines a first phasor of the optical input signal using the first, 

5 second, and third detected powers and a second phasor of the optical 
input signal using the fourth, fifth, and sixth detected powers. Using the 
first and second phasors, the data processor determines the input optical 
signal in the time domain. In effect, the data processing circuitry 
determines from the first, second, and third detected powers a fiirst real 

10 part and a first imaginary part of the input optical signal in a first 
complex reference plane. Similarly, from the fourth, fifth, and sixth 
detected powers, the data processing circuitry determines a second real 
part and a second imaginary part of the input optical signal in a second 
complex reference plane. The first complex reference plane corresponds 

15 to the first polarization state of the optical reference signal, and the 

second complex reference plane corresponds to the second polarization 
state of the optical reference signal. The first real part and the first 
imaginary part correspond to the first phasor, and the second real part 
and the second imaginary part correspond to the second phasor. 

20 Together, the first and second phasors accurately represent all the 

polarization states of the input optical signal and include both the real 
and imaginary signal components in both polarization states. As a 
result, the input signal can be reconstructed very accurately in the time 
domain (more so than when only one phasor is employed). 

25 Although the present invention can be practiced without 

calibration, better results are achieved when cahbration is employed. 
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Typically, optical components will have tolerances that are no better 
than 5% or 10%. Therefore, if the system is constructed from 
commonly available optical components, resulting measurements are 
accurate to this same degree. But higher degrees of accuracy are 

5 required. One method of achieving this would be to buy increasingly 
better optical components, but tiiis approach rapidly becomes cost 
prohibitive. A better way is to accurately calibrate their unwanted 
effects on the measured signals, and then computationally remove them 
from the measurements. In this way, higher accuracy measurements can 

10 be achieved at a much lower cost. In fact, the accuracies tiiat can be 
achieved with this approach may exceed the accuracies that could be 
achieved regardless of optical component cost. 

In a first example caUbration, ampHtude and phase corrections for 
the phasor output by the first detector block may be determined and used 

15 to generate a phasor calibration matrix for the first detector block. As 
the reference signal is swept across a range of different frequencies, 
corresponding detected powers at the first detector block are acquired 
for each frequency. Each detected power has a high frequency 
component. Phase differences are calculated between the detected 

20 powers m the first detector block using the corresponding acquired high 
frequency components. The reference signal may again be swept across 
the range of different wavelengths without an input signal to acquire 
detected powers in the first detector block. Each detected power has a 
low frequency component. Amplitude differences between detectors in 

25 the first detector block are calculated using the corresponding acquired 
low frequency components. These phase and amplitude differences are 
used to generate flae phasor calibration matrix. Preferably, the phasor 
caUbration matrix is converted from an arbitrary reference system into 
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complex plane reference system. A sunilar calibration procedure may 
be applied to the second detector block in the second exam4)le 
embodiment. 

A vector calibration procedure is also desirable in the context of 
5 the second example embodiment described above. Accurate vector 
measurements require that the signal be projected onto two orthogonal 
vectors with the same length. The measurement is accurate only to the 
degree that these conditions are satisfied. If those conditions are not 
satisfied, which they typically are not in a real world system, the system 
10 should be calibrated for the non-ideal aspects of the reference vectors. 
Because optical polarization states within fiber optic networks tend to 
vary over wavelength, (e.g., as the laser source is tuned), keeping the 
two reference states perfectly ortiiogonal is very difficult and expensive. 
If, however, the two states need only be approximately orthogonal, but 
15 repeatable, then the system is much easier to build and only requires a 
calibration procedure to produce the correct vector measurements. 

In the vector calibration procedure, the reference signal is 
generated at multiple different polarizations, (e.g., four), and the 
resulting power phasors determined at the first and second detector 

20 blocks construct a complex vector corresponding to each different 

reference signal polarization. A vector calibration matrix is generated 
using the complex vectors generated for each of the reference signal 
polarizations. The vector calibration matrix is used in normal operation 
to convert subsequentiy detected powers at the first and second detector 

25 blocks into an ortho-normal coordinate system. 
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It is further desirable to calibrate a reference signal generator to 
ensure that the frequency of the generated reference signal matches the 
frequency the generator is set to generate. In the frequency calibration 
procediire, the reference signal is swept across a range of different 
5 frequencies. A portion of the reference signal is passed through two 
different length fiber paths. Light from the two different length paths is 
detected as a function of wavelength. A reference signal frequency 
correction is determined using those detected outputs. 

A final procedure relates to the recovery of the complex spectrum 
10 of the input optical signal. This signal recovery requires the spectrum of 
the signal to remain constant over the course of the measurement, which 
is not difficult if the signal is repetitive with some period that is known 
or can be measured. The processes and calibrations described above 
allow measurement of the time variation of a complex vector with 
15 respect to a known reference field. This raw signal is only a part of the 
overall signal which may have a bandwidth much larger than the 
bandwidth of the detector. By observing the signal over a series of 
bands that completely covers the signal bandwidth, the signal can be 
reconstructed entirely at a much higher effective sampUng rate than 
20 would otherwise be possible. The effective sampling rate is determined 
by the sweep range of the reference laser, and can be on the order of 
1000 GHz. 

A frequency bandwidth response for each detector block is 
determined. From the corresponding frequency response, a time domain 
25 impulse response of each detector block is determined. The impulse 
response is used to create a Green's function that relates the input 
optical signal to be determined as a function of time and the measured 
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signal determined from the detected power as a function of time. The 
Green's function is inverted and then used to convert the measured 
signal into the input optical signal. 

Other features, aspects, and advantages of the present invention 
5 will become apparent from the following detailed description, taken in 
conjunction with the accompanying drawings, illustrating by way of 
example the principles of the invention. Like reference numerals refer 
to like elements throughout. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 FIGURE 1 is a graph illustrating an original modulating signal in 

the time domain; 

FIGURE 2 is a graph illustrating an optical signal modulated with 
the information from Figure 1 in the frequency domain showing just the 
real part of the frequency spectrum; 
15 FIGURE 3 is a graph illustrating a transform of the frequency 

domain signal in Figure 2 into the time domain; 

FIGURE 4 is a graph illustrating the modulated optical signal in 
the frequency domain that includes both the real and imaginary parts of 
the frequency spectrum; 
20 FIGURE 5 is a graph illustrating the original information 

recovered by inverse Fourier transform from the complex frequency 
spectrum shown in Figure 4; 

FIGURE 6 is a depiction of a heterodyne-based optical spectrum 
analyzer (HOS A) in accordance with one example embodiment of the 
25 present invention; 

FIGURE 7 is a depiction of the detector block shown in Figure 6; 
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FIGURE 8 is a flowchart diagram illustrating example 
procedures in accordance with the HQS A illustrated in Figure 6; 

FIGURES 9A-9C illustrate the ou^uts of detector blocks 1 and 2 
shown in Figure 10; 
5 FIGURE 10 is a depiction of an HOSA in accordance with a 

second example embodiment; 

FIGURE 1 1 is a flowchart diagram illustrating example 
procedures for use with the HOSA shown in the second example 
embodiment of Figure 10; 
10 FIGURE 12 is a depiction of a HOSA in accordance with a third 

example embodiment; 

FIGURE 13 is a function block diagram illustrating example 
procedures for phasor calibration of detector blocks 1 and 2; 
FIGURE 14 is a flowchart diagram illustrating example 
15 procedures for vector calibration of the output of detector block one 
projected onto one polarization plane of the reference signal and the 
ouQ)ut of detector block 2 projected onto another polarization of the 
reference signal; 

FIGURE 15 is a flowchart diagram illustrating example 
20 procedures for frequency calibration of the reference laser; and 
FIGURE 16 is a flowchart diagram illustrating example 
procedures for complex spectrum recovery of the input optical spectrum. 

DETAILED DESCRIPTION 

The following description, for purposes of explanation not 
25 limitation, sets forth specific details, such as particular components, 

electronic circuitry, techniques, etc. in order to provide an understanding 
of the present invention. But it will be apparent to one skilled in the art 
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that the present invention may be practiced in other embodiments that 
depart from these specific details. In other instances, detailed 
descriptions of well-known methods, devices, and techniques, etc. are 
omitted so as not to obscure the invention with unnecessary detail. 

5 Individual function blocks are shown in the Figures. Those skilled in 
the art will appreciate that functions may be implemented using discrete 
components or multi-function hardware. Processing functions may be 
implemented using a suitably progranmied microprocessor or general 
purpose computer, using an application specific integrated circuit 

10 (ASIC), and/or using one or more digital signal processors (DSPs). 

A first example embodiment is described in conjunction with the 
heterodyne-based optical signal analyzer (HOSA) 10 depicted in Fig. 6. 
The optical input signal to be analyzed Es = pse'®^^^ (hereafter sometimes 
referred to simply as Eg) is provided to an optical coupler, which in this 

15 embodiment may be a 3 x 3 coupler 12. Any coupler may be employed, 
and one non-limiting example is Gould part no. 23-40355-33-01201 
manufactured by Gould Fiber Optics Division of Gould Electronics of 
Baltimore, Maryland. A laser reference signal Er = ae'^^^^ (hereafter 
sometimes referred to simply as Er) is provided to another input terminal 

20 of the 3x3 coupler 12. The two signals Es and Er are mixed in die 

coupler 12 which provides three interference output signals Dl, D2, and 
D3. The interference signals Dl, D2, and D3 are provided to an optical 
detector block 16 which detects the corresponding power level of tzch 
of the detected signals Dl, D2 and D3 and generates three 

25 corresponding digital power signals PI, P2, and P3. The detector block 
16 may be any suitable detector block, and one non-limiting example is 
a Thorlabs PDA 400 optical detector manufactured by Thorlabs of 
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Newton, New Jersey. The digital signals PI, P2 and P3 are processed 
by a data processor 20 to calculate the original optical signal to be 
analyzed in the time domain Es(t). 

Figure 7 shows an example of the detector block 16 which 
5 includes three power detectors 40a, 40b, and 40c. Each power detector 
40a, 40b, and 40c includes a photodetector 42, coupled to an ampUfier 
44, coupled to a low-pass filter 46. The filtered output is converted into 
a digital format by digital-to-analog conversion means 48, and the 
digital signal is stored in the buffer 50 before being processed by the 
10 data processor 20. 

Figure 8 is a flowchart diagram illustrating the basic steps for 
determining the input optical signal in the time domain. An optical 
reference signal and the input optical signal to be determined are mixed 
to generate multiple mixed signals (step S 1). Multiple power signals are 
15 detected from the mixed signals (step S2). The input optical signal is 
determined in the time domain from the multiple detected power signals 
(step S3). The mathematical procedures implemented by the data 
processor 20 to determine the input optical signal from the multiple 
detected power will now be described. 

20 Given two optical signals incident on a 3x3 detector, where the 

two fields are given by, 

E^^ae'^ (1) 

and 

E,=p,e^' (2) 

25 

the three detected powers will be, 
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Pi=a^ + p^ + 2ap,cos\fi-e+<Pi] (3) 

P,=a^ + p; + 2ap, cos[fi -d+q}^] (4) 

Pi=a^ + p^ + 2ap, cosCS - e + 9)3 ] (5) 

5 where a is the amphtude of 4, is the phase of p is the amplitude 
of Es, 6 is the phase of Eg, and <Pi , <p2 . and are phase shifts detected 
in the 3x3 coupler 12. These phase shifts can be derived using 

conservation of energy and are found to be — for a 3x3 coupler. In 

other words, (p^, (p^, and (p^ are approximately and respectively 0°, 120°, 
10 and 240". 

The signal that we wish to recover is 

E, = p,e'^-'^ = p, cosOS -e)+p, sm{P -9). 

If we can know d , the phase of the reference signal Er, the original 
signal's electric field can be reconstructed as follows. 
15 We begin by subtracting P3 from and P,, 

Pi - P3 = 2ap, cos[/3 -9 + ^J- 2ap, cos[fi - 0 + flJj] 
(6) 

P, - P3 = 2ap^ cos[/3 -9 + (p2]- 2ap, cos\fi -0 + (Pj]. 
20 (7) 

We then divide through by a , (measured previously so it is known), and 
expand the cosines of angle sums into products of sines and cosines. 
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^» = cos()3 -0)cos((p,)- p, sinOS -0)sm(<p,)-p, cos(j8 -0)cos(<P3) 
2a 

+ p,smO3-0)sin(<P3) (8) 



2a 

+ p,sinO3-0)sin(<P3) (9) 



^— ^= p, cosOS -0)cos(ipj)- p, sm(p-d)sm{<P2)-p, cos(p -e)cos(<Pi) 
2a 



5 Collecting terms produces 



10 



15 



^5—5. = p cosOS - 0)[cos(«p, )- cos((jaj )]- p, sin(fi - 0)[sin(p, )- sinCijJj)] 
2a 

(10) 

= cosOS - dfcosM- cos(<P3)]- p, sm(j8 - 0Xsin(<P2)- sin(% )] 

2a 

(11) 

We can rewrite the two equations above in vector notation 



p —P^ 
2a 



2a 



[cos((Pi)- cos(<P3)] [sin(«P3)-sin(<p,)]Tp, cos(^ 
[003(^3)- cos(p3)] [sin(<P3)- sin((Pj)]jLp, sin(fi -e)\ 



(12) 



The vector 



p,cosO3-0)' 
_p,sin(j3-a) 



(13) 



is the vector in the complex plane of the signal Eg to be reconstructed. 
20 We can invert the matrix and multiply both sides of the equation to find 
this vector. 
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■ [cos(<j!)i )- cos(<P3 )] [sin {(p^ )- sm{(p^ )] 
(14) 



M -'a 

2a 



2a 



p,cos(i3-6>)" 
p,sinO8-0)_ 



For illustrative purposes, assume that the 3x3 coupler is ideal, and, 

<P,=0 (15) 

(16) 

(17) 



27C 

f2 3 



2w 



10 We can then evaluate the sines and cosines exactly to get, 





-0 


2 

LV 





2 



2a 



2a 



(18) 

IS and with simplifications. 



2 




2 


2 


0 




(19) 





-1-1 




p^cos{P-9) 



20 Inverting the matrix gives. 



2 

0 - 



(20) 



73 



P " P 

2a 
P - P 

2a 



p,sin(j8~0) 
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10 



and multiplying through gives. 



3a 6a 



S[ 2a 



(21) 



As a result, the desired electric field Ej can be expressed as follows 
E, = P,e'^P-'^ = p, cosOS -d)+ip, sinOS - 0) = 



3a 6a J VSl 2 



(22) 



p — P P - P^ 
3a 6a 
(23) 



1 fp,-p. 



V3 



2a 



The data processor uses the known reference signal amplitude a and the 
15 known detected powers P1-P3 to determine the input optical signal Es. 

The first example embodiment assumes that the optical signal to 
be determined Eg has only one polarization. But as mentioned in the 
background, optical signals have more than one polarization. If only 
one polarization state is detected, as described so far, then there is a 

20 possibility that the unknown signal could be polarized orthogonally to 
the detected state. In that situation, no signal would be detected. A 
more likely case is that the signal will he somewhere between perfect 
ahgnment, (i.e., in the same polarization state), and perfect 
orthogonality, (i.e., zero detection). The signal will be detected, but not 

25 its absolute amplitude (power) because some percentage of the signal 
will be in the orthogonal state and will go undetected. A further 
complication occurs if the signal changes polarization as a function of 
time or frequency. In this case, it will be impossible to detect this 
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polarization change. Instead, the polarization change will be wrongly 
interpreted as a change in amplitude or phase. 

To address these shortcomings and account for the two 
polarizations of the input Ught signal, the input optical signal is treated 
5 in the second example embodiment as a vector quantity with two 
orthogonal vector components called first and second phasors. These 
two phasors fully describe the optical signal to be determined. Each of 
the first and second phasors is a complex signal having a real part and an 
imaginary part. 

10 This concept is represented graphically in Figs. 9A-9C. Figure 

9A shows vector Eg which represents the complex optical signal that to 
be reconstruct in the time domain for analysis. The vector Es is 
measured with respect to two orthogonal polarization planes pl(t) and 
p2(t) of the reference signal Ef. The first polarization plane Er,pi(t) 

15 corresponds to the vertical axis in Fig. 9A. The second polarization 
plane Ef,p2(t) corresponds to the horizontal axis in Fig. 9A. The vertical 
component of Eg corresponds to the dot product of vector Es and the 
vertical polarization plane E^^pi. In other words, the first phasor is the 
projection of vector Es onto the vertical polarization plane Er^pi. 

20 Similarly, the second phasor is the Eg vector corresponds to the dot 
product Es-Er,p2 and is the projection of Eg onto the horizontal 
polarization plane. 

The vertical component of vector Es is the first phasor Pi 
illustrated in Fig, 9B. Phasor PI corresponds to the polarization plane 
25 Er,pi and has both an amplitude Api(t) and an angle of direction 0pi(t). 
Phasor PI is a complex number and has both a real component 
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corresponding to the horizontal axis and an imaginary component 
corresponding to the vertical axis. Similarly, Fig, 9C illustrates the 
second phasor P2 which corresponds to the horizontal component of the 
vector Es in the polarization plane Er,p2. It also is a complex number 
5 having an amplitude Ap2(t) and a phase 6pi(t) and real and imaginary 
components. 

In order to account for the two polarization plane components of 
the vector Eg, two detector blocks are employed in a second embodiment 
of an HOSA shown in Fig. 10. The detector block 16 corresponds to 

10 phasor PI shown in Fig. 9B, and the detector block 30 corresponds to 
phasor P2 shown in Fig. 9C. The two reference signal polarization 
planes Er^pi and Er,p2 are generated by splitting the laser reference signal 
Er in a two-by-two splitter 24. One output of the two-by-two splitter 24 
corresponds to the reference signal in an Er,p2 polarization which is 

15 forwarded to a three-by-three coupler 28. The reference signal in the 
other output from the two-by-two splitter 24 is provided to a 
polarization controller 26 which changes the polarization of the 
reference signal from Er,p2 to Er,pi which is approximately orthogonal to 
polarization Er,p2. One non-limiting example of a commercially 

20 available polarization controller is the FPC031 from Thorlabs. The 
reference signal Er,pi is provided as an input to the first three-by-three 
coupler 14 as described above. The second three-by-three coupler 28 
provides three outputs of the mixed signals Eg and Er at the polarization 
Er,p2 as inputs D4, D5 and D6 to detector block 2 (30). Detector block 2 

25 detects a power signal P4, P5, and P6 corresponding to the detected 

nodxed signals D4, D5 and D6 and provides those detected power signals 
to a data processor 32. 
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Data processor 18 processes the three power levels to determine 
the real and imaginary components of the phasor Pi, and data processor 
32 determines the real and imaginary components of phasor P2. Recall 
that the phasor PI is the projection of the input signal to be determined 

5 Es onto the reference signal Er in polarization PI , and phasor P2 

corresponds to the projection of the signal Eg on the reference signal in 
polarization P2. The complex numbers corresponding to phasors PI and 
P2 are provided to data processor 20 which determines the input signal 
to be analyzed vector Es shown in vector form in Fig. 9A. Data 

10 processor 18 performs the calculations in equation (24) set forth below; 
data processor 32 performs the calculations in equation (25) set forth 
below, and data processor 20 performs the calculations in equation (26) 
set forth below. 

Fig. 11 is a flowchart diagram illustrating example procedures in 
15 accordance with the second example embodiment shown in Fig. 10. The 
optical reference signal is first split into fnst and second reference 
signals in splitter 24 (step SI). The input signal is also split into first 
and second input signals in splitter 12 (step S2). The polarization of the 
first reference signal from the splitter 24 is changed in polarization 
20 controller 26 (step S3). The first reference signal and the first input 
signal are mixed to generate first, second, and third mixed signals 
corresponding to Dl, D2, and D3 (step 84). The second reference signal 
and the second input signal are mixed to generate fourth, fifth, and sixth 
mixed signals corresponding to D4, D5, and D6. respectively, (step S5). 
25 Detector blocks 1 and 2 detect first through sixth power signals from the 
first through mixed signals (step S6). 
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Data processor 18 determines from the first through third 
detected power signals a first real part and a first imaginary part of the 
input signal Es projected into a first complex polarization plane PI (step 
S7). The data processor 32 determines from the fourth through sixth 
power signals a second real part and a second imaginary part of the input 
signal Es in a second complex polarization plane P2 (step S8). Data 
processor 20 calculates the input signal Es in the time domain using the 
first and second real and imaginary parts (step S9). 



10 



15 



The mathematical processing performed by data processor 20 is ' 
now described. The same procedures may be performed on detector 
block 2 as described earlier for detector block 1 to obtain the electric 
field parallel to the other reference signal polarization. We now have 
two projections of the electric field Es, and thus the electric field vector 
E, , in some arbitrary x, y plane space: 



sy 



3a, 



6a. 



— I 



V3 



6a^ 



P2-P3 
2a 



P5-P. 

2a 



(24) 
(25) 



20 



The electric field vector may be converted to an ortho-normal space by 
multiplying this vector E^E^ by some transforming matrix, T , 



E,=T 



(26) 



3a„ 6a„ 



3a 



6a, 



— I 



— I 



V3 



P. -P. 



P -P 

2a„ 
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The notation may be simplified by letting the complex phasor 
measurements be simple dot products and by not doing any 
normalization in the first step. We start again with. 



(27) 
(28) 



V31 2 J 

where E^^^ , and , ^ reference fields at detector blocks one and 
two, respectively. Simplifying these equations further, 

10 -^.1 =7(2^1+^2 -^P3)-i:;^iP2-P^) (29) 

(30) 



2^3 
1 

2V3 ' 



allows rewriting the vector as: 



E^=T 



{2P, + P,-3P,)-i^{P,-P,) 
_{2P, + P,-3P,yi^{P,-P,\ 



(31) 



Equation (3 1) is a full vector representation of the electric field so 
that power or phase measurements are not polarization-dependent. Also, 
polarization dependent signal characteristics (such as the effects of 
polarization mode dispersion) are measurable. 

A third example HOSA embodiment is depicted in Fig. 12. 
Various calibration procedures will be described in the context of the 
HOSA in Fig. 12. However, some of the calibration procedures may be 
employed with other embodiments. A caUbration mode switch 62 
coupled to a second polarization controller 64 are inserted between the 
optical input signal to be analyzed and the two-by-two sphtter 12. When 
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the mode switch is set to "1," the input signal Es is delivered to the 
polarization controller 64. When the mode switch is set to "0," the 
reference signal Er is sent to the polarization controller 64. Also 
included is structure for calibrating the tunable laser 22 which includes a 
5 three-by-three coupler 66 receiving the reference signal from the splitter 
24, receiving inputs from two Faraday Rotating Mirrors (FRMs) 68 and 
70, and generating outputs detected by a third detector block 72. 

A first calibration procedure is the adjustment of polarization 
controllers 62 and 64. With the mode switch 62 in the "0" position and 

10 the reference signal 22 (the output of the tunable laser) in a continuous 
sweeping mode being incremented through a particular frequency range, 
such as 1550 nm to 1552 nm, the polarization controller 64 is adjusted to 
maximize the power levels detected by detector block 30. When this is 
accompUshed, the power level at D5 in detector block 30 is minimized 

15 by adjusting polarization controller 26. At that point, polarization 

controller 64 is adjusted so that the power levels at P2 in detector block 
16 and P5 in detector block 30 are approximately the same (within plus 
or minus 10%). This calibration procedure ensures that the two detected 
polarization states are not the same and are approximately orthogonal, 

20 which leads to better quality measurements. 

A next calibration relates to a phasor calibration of detector 
blocks 1 and 2 (16 and 30). For each three-by-three coupler 14 and 28, 
the coupling coefficient from any coupler input to any coupler output is 
equal, and the power detected on each of the outputs should be separated 
25 in phase by exactly 120°. In reality, the coupling coefficients between 
coupler outputs are not the same. Nor is the phase difference between 
each of the three detector outputs in each detector block exacdy 120''. 
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The phasor calibration procedure compensates for these coupling 
coefficient and phase shift differences. 

The phasor calibration begins with a series of measurements 
aimed at assigning a complex number z^-EJE^, which will represent 
5 the measurement, to detector block 1 . After the two polarization 

controllers are set properly, the mode switch is set to the "0" position. 
The optical power signal P measured at the n^ detector is 
= a'r^ + ik-aX + 2&ni,r„a=^ cos[^(f)- -r)- 



10 where n=l, 2, 3 denoting tiie three detectors Dl, D2, and D3 in Detector 
Block 16. The amplitude of the reference signal is a . The field loss 
from die laser 22 through tiie splitter 24, the PC 26, and the coupler 14 
to the detector "n" is r„ . Since the cosine term (i.e., the AC term) 
oscillates rapidly as the reference laser is tuned, this term can be 

15 extracted from the equation via the following sequence: 1) Fourier 
transform to the time domain, 2) window in the time domain to 
separate the DC terms from the cosine (AC) term, and 3) transf oral back 
to the frequency domain. This sequence 1-3 leaves us the two 
quantities: 



(32) 



20 



(33) 



and 



iPn\c = 2^aX^n cos[0(r) - 0(f - tr) - 



(34) 



25 



The pa terms can be found from Eqs. (35)-(37) 
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(35), (36), (37) 



where, since we are interested only in relative phase, we have chosen to 
reference and to and set A = 0 • With the mode switch set to 
the "1" position and with no signal present, we measure at each of 
5 the n detectors. Then with the measurement signal present, we measure 
at each of the n detectors. These power signals will be of the form 

(38) 

and 

10 P„'" = + p^m^ + 2apm„r„ cos[0(O -0(0 - jSJ. 
Using the following definitions: 



pwi pc pm pc 

p -._£2__£2 M M 



A = 



B = 



pm pc p/n pc 

p — _£3__13 M 

H[Pzli0C-P2 



P' 



■sin[i3J-, 



7COs[A], 



(39) 



(40)-(45) 



sin[j8i]. 



C = 



{p'-\oc-pr "iiPXc-Pi 



P' 



sinEAl, 



we get the linear equation 



wo 2004/005974 



PCT/US2003/021337 



27 



= 2 



C D 



ka 
ka 



cos(0 - ^) 



sin(0-^) 



(46) 



the solution of which is 



-^cos(0-0) 
ka 



^ ka 



sin(0-0) 



2 C D 



-1 , 



(47) 



In order that this solution exists, the matrix must be invertible. 
Using Eq. (47), along with the expression for ^„ given in Eqs. (35)-(37), 

we can assign z, = -^[cos(0 — ^) + i'sin(0 —0)] ,to detector block one. 

ka 

10 Following the same procedure, we can 

assign = — [cos(9 - ^) + i sin(0 - 0)] to detector block two. The 

' ka 

calibration of the two detector blocks can be done simultaneously. 
Finally, the difference in the relative phase of the complex detectors as a 
function of the laser frequency must be included in a constant phase 
15 term that will be included in the matrix. The phase term will be a linear 
function of frequency and results from mismatches in the interferometer 
lengths. We then end up with an imcalibrated vector measurement of, 



20 



1 



ka 
(48) 



^ ;,y(«,('W('>r.) 



ka 



These calculations are made in data processors 18 eind 35, and the 
two components are assembled into a vector in data processor 20. The 
vector calibration matrix described below converts Eraw to an accurate 
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representation of the electric field in an ortho-normal basis set. We can 
then use this measurement in the time-domain reconstruction or other 
calculations such as simple spectral power or polarization variation as a 
function of wavelength. 

5 The phasor calibration of detector blocks 1 and 2 will now be 

described in more general terms in conjunction with the flowchart 
diagram depicted in Fig. 13. The mode switch 62 is set to zero to 
calibrate the reference signal (step SI). The reference laser is swept 
over a desired frequency range, and at each frequency in the range, 

10 powers having both DC and AC components are detected at each of the 
detector blocks 1 and 2 as a function of wavelength (step S2). Using the 
detected AC power at each detector block, relative phase differences are 
calculated between each of the three detectors in each detector block 
(S3). 



15 ^ The mode switch 62 is set to "1" with no input signal, i.e., there 

is no signal input into the polarization controller 64 (step S4). The 
reference laser is swept again over the desired frequency range, and the 
resulting DC power for each of the detector blocks 1 and 2 is determined 
for each frequency scan (step S5). Relative amplitude differences 

20 between each of the three detectors in each data block are determined, 
i.e., Pi-Ps and P2-P3 for detector block 1, and P4-P6 and P5-P6 for 
detector block 2 (step S6). A phasor calibration matrix is generated for 
each of the detector blocks using the DC powers and relative phase 
differences for each detector block as detected in steps S3 and S6 (step 

25 S7). The phasor calibration matrix for each detector block is converted 
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from an arbitrary x, y reference system into an accurate phasor 
calibration matrix in the complex plane (step S8). 

The phasor calibration allows for the correct assignment of a 
complex number to the electric field incident on the detector block— 
5 through the determination of the constants A, B, C, and D— that 

represents the measured signal to detector blocks one and two. Without 
this calibration, accuracies better than 10% would be very difficult to 
achieve. With this calibration, accuracies on the order of 1% should be 
achievable, 

10 Another desirable calibration is referred to as vector calibration. 

Vector calibration is designed to make the system insensitive to the 
polarization of the measured signal. The vector calibration has the 
further benefit of providing a calibrated measurement of the polarization 
of the incoming light so that effects like polarization mode dispersion 

15 (PMD) can be measured and observed. 

The vector calibration begins by setting the mode-switch to the 
"0" position and taking measurements at detector blocks one and two for 
four distinct settings of the polarization controller PCa- Using the 
results of the vector calibration, a complex number can be assigned to 
20 bodi detector block one and two for each of the four measurements. At 
detector block one, we get Zn, Zi2, Zo, z^, and at detector block two we 
get Z21, Z22, Z23, Z24. With these complex numbers, we form the 
following vectors: 
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V3 = 















( 7 ^ 




'Zu^ 






[Z^j 



(49)-(52) 



With these definitions, the following matrix can be formed 



'P 8^ 



(53) 



where [3c y]is a matrix with colmnns formed by the elements of the 
vectors x and y . Using the following set of definitions 

W = \pf + \qf-U 
X=\gf+\hf-l, 
Y=\pq\, 

z=\ghl 



r = Arg 



WZ + XY .WZ-XY 



-hi- 



cosy 



siny 



sin 0 = 



qh] 

-W 



2Y cos {Arg[p/q]+^y 
cos © = Vl — sin^Q, 



(54)-(62) 



the vector-calibration matrix is given by 

[V, vj]-'. (63) 



M = 



^1 e'^sine"^ 



0 COS0 



With the switch in the "1" position, a measurement vector is 
formed by making a measurement with the signal present and applying 
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the scalar calibration. With the mode switch 62 in position 0, which 
directs a portion of the tunable laser signal to polarization controller 64 
and the subsequent elements, the tunable laser is swept, and a 
measurement of the signal generated is made. It is presumed that the 
complex detector blocks 14-18 and 28-32 have been previously 
calibrated and produce complex numbers representative of tiie 
projection of the electric field that went through the measurement patii 
62, 64, 12, etc. onto the reference electric fields entering each complex 
detector block, Er,pi and Er.p2. These two complex field measurements 
firom the two complex detector blocks are the entries in a vector that 
characterizes the electric field of the light. This vector field is in an non- 
orthonormal basis set and must be converted into an orthonormal basis 
set before accurate measurements can be obtained. 

This leaves us with two complex numbers Zim and Z2m, where, 
again, the 1 and 2 refer to data taken at detector blocks one and two, and 
the "m" refers to measurement (as opposed to caUbration). The 
polarization calibrated measurement vector is then given by 



This calculation occurs in processor 20, and the accurate vector 
measurement produced is used to calculate the final measurements of 
the HOSA. Among these measurements are the time-domain response 
of the system, the amplitude spectrum, and the polarization state as a 
function of time or frequency. 




A general overview of the vector calibration of detector block 1 
(vector projected onto the first polarization plane of the reference signal) 
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and the data block 2 (vector projected on the second polarization plane 
of the reference signal) is now described in conjunction with the 
flowchart diagram in Fig. 14. Mode switch 62 is set to "0" so that the 
reference signal is provided to the polarization controller 64 (step SI). 
5 The polarization controller 64 is regulated to four distinct positions by 
processor 20 that correspond to four different polarization angles (step 
S2). For each of the four polarizations, the complex numbers detected 
by the phasor-calibrated detector blocks 1 and 2 are recorded (step S3). 
A vector calibration matrix is constructed from the four recorded 
10 complex numbers to convert subsequentiy detected complex numbers 
from the first and second detector blocks onto an ortho-normal plane 
(step S4). Without this calibration, accuracies better than 10% would be 
very difficult to achieve. With this calibration, accuracies on the order 
of 1% should be achievable. 

A further desirable calibration of the HOS A is a frequency 
calibration of the reference laser signal. The measured complex vector 
representing the input optical signal Es is an instantaneous measurement 
in time, and thus, contains littie information about the spectrum of input 
optical signal Ej. That spectrum is reconstructed by sweeping the 
reference signal across the frequency spectrum of the input optical 
signal Es. Thus, it is important that the reference signal generator/lasCT 
be very accurately calibrated. 

Reference is now made to the frequency calibration of the 
reference laser E^ flowchart shown in Fig. 15. The reference laser is 
25 swept over the frequency range of interest corresponding to the input 
optical signal (step Si). A portion of the reference signeil is passed by 
way of the three-by-three coupler 66 to each of the Faraday Rotating 



15 



20 
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Mirrors 68 and 70 through different length fiber paths (step S2). The 
Faraday Rotating Mirror 70 has a longer fiber path than the Faraday 
Rotating Mirror 70. The reference signal is reflected back by its 
respective Faraday Rotating Mirror to two inputs of the three-by-three 

5 coupler 66. Because of the different fiber lengths, there is a delay 

between the two reference signals received at the coupler 66 from (step 
S3). Because of the different fiber lengths, there is a delay when the 
reference signal is received the coupler 66. The detected power levels 
on detector inputs D7 and Dg of detector block 3 (72) are recorded (step 

10 S4)- This occurs for each swept reference laser frequency. The average 
and peak-to-peak amplitude variations at detectors 7 and 8 are 
determined (step S5) as well as the phase difference between detectors 7 
and 8 for each swept frequency (step S6). The amplitude and phase 
variations are used to determine the actual frequency output of the 

15 reference laser. 

The measured phase of the signal Eg is the difference between the 
phase of the reference Er and the phase of the input signal Eg. The phase 
of the reference Er can be obtained from the reference wavelength 
monitor shown at the bottom of Figure 12 and consists of a 3x3 coupler 

20 66, two Faraday Rotator Mirrors 68 and 70, and detector block 3 (72). 
The two signals on detectors D7 and D8 of detector block 3 give the sine 
and cosine of a phase that is proportional to the frequency of the 
reference Ef. The signal on detector D9 is proportional to the power 
offset on the sine and cosine signals. This offset power can be 

25 subtracted off either electronically or optically. The remaining sine and 
cosine terms contain a static relative phase error between the signals. 
This phase error can be corrected for in the following way. 
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There are two power readings, 
= + by cos (on + Si) 



(65) 
(66) 



where a, and are the offsets of the fringes on the respective 
detectors, and and by are the respective fringe amplitudes. The phase 
shift Q. is the phase difference between the two channels and will be 
approximately — , where the coupler is sn NxN coupler. We are 
looking for the angle, cor , which can be found to be, 



(or = Arg 



b. 



J J 



Once the phase signal proportional to optical frequency has been 
obtained, it can be integrated to recover the original phase of the 
reference source. Since the complex numbers Px, Py have a phase that is 
relative to the reference Er phase, division by a complex number having 
the phase of the reference Er a number having only the phase as the 
signal of interest: 



20 ^{t)^^\o){t')rdt\ 



(68) 



This is the phase of the reference signal, which appears as j3 
in the discussions of the detector cahbrations above. Given knowledge 
of the reference phase from equation 67 and a measurement of the signal 
25 phase difference from equations (31) or (48), we can calculate the phase 
of the signal field, Because of the limited bandwidth of detectors 
D1-D6, the signal is measured over a limited (e.g., -lOMHz) bandwidth 
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at any given tune. The naeasured signal as a function time is thus only a 
portion of the overall signal of interest. By sweeping the reference laser 
over the full range of the signal, all of the signal can be observed. Since 
the bandwidth of the detectors and the precise frequency of the reference 
5 laser are known, the full time-domain signal may be reconstructed 
analogous to assembling a panoramic view by joining together a set 
photographs that slightly overlap. 

In this case, however, the signal is varying with time, and the 
snap shots (continuing with the panorama analogy) must be 

10 synchronized. This synchronization is achieved by sampling the signal 
at an integer multiple of the signal repetition rate. Although the 
assembling of snap shots is a good heuristic picture of the process, the 
optimal practice differs substantially. We have the equivalent of a very 
precise measurement of the camera angle (the reference phase) and 

15 fading on the picture edges (the bandwidth limitation on the detectors). 
In order to properly recover the full panoramic view, we must blend 
overlapping pictures, weighting each appropriately by the degree of 
fading on each at any particular pixel. 

What follows below is a mathematical description of a process 
20 for assembling the overall time signal. The mathematical nature of the 
description means that it is easy to translate the analysis into a computer 
algorithm. Thus the description below is easy for the computer to 
understand, and the description above is easy for the reader to 
understand. In this regard, the analysis below should be regarded as a 
25 particular example solution. Many other solutions may be found. 
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Having a particular solution at hand merely illustrates that there is such 
a solution, and it can be implemented on a computer. 

We begin with our detected and calibrated complex signal, 
5 zM-[zA^)e"^'^]®bit) (69) 

here z„(t) is the measured complex signal at the detector block(s) and 
z^(^)is the original input signal that we wish to recover, 0(r)is the phase 
of the reference of which we have a measurement, and b(t)is the impulse 
10 response of the band-limited detector or block 1 (and block 2?) and ® is 
the convolution operator. We can then take a Fourier transform to get, 

J z„ (r )^'-rfr = B{q))] z, (T)e'^^^^^«rfr. (70) 
Invert the transform gives 

15 

]]z^ {r)e'^dre'^^d(o =] B{co)] z, {r)e''^''^^'^dre''^'d(0. (71) 

Collecting terms under a double integral, 

20 z^ (0 = J J (r )5 Wco, (72) 

and extract the unknown signal function from the integral over 
frequency yields 

25 z^ (0 = J 2. (^) J ^ (cy )e^^('^^-'"('-^W£t>jT. (73) 

We can the define a new function of two time variables, forming the 
Greens function that relates the measured signal to the input signal 
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b'{t,x) = / 5 {(o)e''^'V'"^'-'^d(o. (74) 

Note that this function is comprised entirely of known or measurable 
5 quantities. Substituting back in, we then get 

zM=]z,{r)b'{t,x)dT. (75) 

So, the measvurement is simply the result of a linear operator on the 
10 signal, z, having a kernel given by b'if,x) . If the operator is 

invertible, then we will be able to reconstruct the signal. If we discretize 
the process, Eq. (75) becomes 



15 



Zk = z„ {kt) = {nr)b\kt,m)^T. (76) 
Assigning new names to the discretized functions gives 

^*=i;«AAT. (77) 

/i=0 

20 Rewriting in linear algebra notation: 

z = te. (78) 

If the calculated matrix, b , is invertible, then the discrete signal, a , can 
25 be recovered by, 

(79) 

r? = a (80) 
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In general, ^is a very large matrix (--1 million by Imillion), however, in 
most cases it will be extremely sparse having only 10 to 100 elements 
centered on the diagonal. This should make the numerical problem 
tractable. 

5 The recovered vector, 5 , is an accurate reconstruction of the 

electric field as a function of time, and from this signal, any 
characteristic of the optical spectrum of the original input optical signal 
can be recovered. 

A brief overview of the complex spectrum recovery is now 
10 described in conjunction with the flowchart depicted on Fig. 16. The 
electrical frequency response for the bandwidth of the two detector 
blocks 1 and 2 is determined (step SI). The time domain impulse 
response of each detector block is then calculated from its frequency 
response using the inverse Fourier transform (step S2). The impulse 
15 response to create a Greens function that relates the actual input signal 
as a function of time to the measured input signal as a function of time 
(step S3). The Greens function is then inverted to find a function that 
will convert the measured input signal generated by the detector blocks 
1 and 2 to the actual input signal (step S4). In the context of the 
20 aboveanalogy, this step is the blending of the large set of snap shots that 
produces the overall view. It is however described in mathematical 
terms than one skilled in the art of algorithm creation can imderstand 
and directiy put to use. 



25 



While the invention has been described in connection with what 
is presentiy considered to be the most practical and preferred 
embodiments, it is to be understood that the invention is not to be 
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limited to the disclosed embodiments, but on the contrary, is intended to 
cover various modifications and equivalent arrangements included 
within the scope of the appended claims. 



